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TECHNICAL PUBLICATION 


STATISTICAL EVALUATION AND IMPROVEMENT OF METHODS 
FOR COMBINING RANDOM AND HARMONIC LOADS 


1. INTRODUCTION 


Many structural components in dynamic environments experience both sinusoidal and random 
loadings at the same time. One example of such a structure is a component of a rocket engine. There are 
several sources that produce the random component of the excitation, including the gas generator and main 
combustion chamber. These excitations are characterized by a power spectral density (PSD) spectrum of 
accelerations. These accelerations are used in a base excitation random response analysis, which results in 
a PSD response spectrum at any desired location, as shown by the example of the gas generator axial force 
response for NASA Marshall Space Flight Center’s MC-1 engine in figure l. 1 Although it can be in any 
desired quantity, the mean of this spectrum, defined simply as a load, is zero, and since the process can be 
assumed to be Gaussian, 2 the root mean square (rms) of the response is defined to be equal to one standard 
deviation, a. Simultaneously, the turbomachinery in the engine generates large harmonic loads due to the 
unavoidable unbalance in their rotors. The frequency of the harmonic excitations are at the shaft speed of 
the turbines and their multiples. The excitation from this sinusoidal loading, which can also be determined 
by previous acceleration measurements, is applied independently in a frequency response analysis, which 
generates the amplitude of a sinusoidal response for any desired location for a specified range ol harmonic 
excitation frequencies. 



Figure 1 . Typical results from random response analysis. 



2. LOADS COMBINATION EQUATIONS 


In general, the structural dynamic response analyses for the random and harmonic cases are 
performed separately, and the results are then combined to obtain a design load. These loads are used to 
determine a level of alternating stress in the component for design to ultimate or yield strength and for 
high-cycle fatigue (HCF). The calculation of this load has been a subject of some dispute in industry. 
Although the development of methods for combining different types of random loads has been a subject of 
extensive research, 3 detailed study of the use of similar equations for the combination of random and 
harmonically generated loads has not been performed. Perhaps because of this lack of examination, an 
industry standard loads combination equation has not been chosen. There are, therefore, several different 
methodologies that have been applied by various users. They are all based on the method used for obtain- 
ing design loads for purely random environments. Although there have been studies suggesting that the 
extreme-value distribution should be used for limit load design, 4 a widely accepted practice is to use 
3 times the rms of the random load, < 7 ran , for the design load where the distribution is assumed to be 
Gaussian. This value will exceed the response 99.865 percent of the time. It has also been suggested that a 
value less than 99.865 percent is acceptable for HCF design, but this Technical Publication (TP) will 
assume that both the Gaussian distribution and the 99.865-percent excedance level are valid. 

The most frequently used method for combining the two dynamic load components, referred 
to here as the standard method, simply uses the sum of the amplitude of the sine load A and 3xa : 5 

sin ran* 

design load = A sin + 3cr ran . ( 1 ) 

This technique assumes the sine peak value always occurs simultaneously with the peak random value. 
Another frequently applied method is to multiply 3 by the square root of the sums of mean square values of 
both the random and the sines, called the 3xssMS method: 6 

design load = 3^(a sin ) 2 +(<T ran ) 2 , (2) 

where 0" s | n , the rms of the sine wave, is equal to the sine amplitude divided by -^2 • This method treats the 
sine as if it had a normal probability distribution, which has a much larger range than a sine function. The 
3xssMS approach is similar to equations written to combine loads on launch vehicles, which do not have a 
substantial harmonic component but do have random loads from several sources. 7 A third approach, 
proposed by Steinberg, 8 is to take the root sum square of the peak values (called the peak method): 


design load = 




4 '(3cr r an)~ 


(3) 


2 


Although it is stated in reference 8 that the standard method is too conservative and that it is convenient to 
use the suggested peak method, no mathematically rigorous justification is presented for using it for either 
ultimate/yield or HCF design. None of these closed-form methods produce a design load that can be con- 
sistently associated with a particular probability level. 

Although the methodology discussed in this TP is relatively straightforward, it has not been gener- 
ally applied in industry. In a paper by DebChaudhury et al., advanced probabilistic techniques were used to 
evaluate the stochastic nature of the parameters of the loading and the structure for the Space Shuttle main 
engine (SSME), but the sine and random loads were combined using the simple closed-form 3xssMS 
method. 2 A similar study by Newlin et al., in which stochastic techniques were used to evaluate HCF 
problems in the SSME, generated a simulated random time history and added it to a harmonic time 
history. 9 The combined loading was then used for each time step. This methodology is similar to that 
presented in this TP but is too complex to use in a routine manner. Schock 10 was apparently able to analyti- 
cally calculate the joint cumulative density function for (7=1 and was able to generate charts for different 
cummulative distribution function (CDF) values for this level of random a. Ferebee 1 1 added to this work; 
however, the results do not appear to be consistent with the work presented below, and the methodology for 
performing the integration is not discussed. 
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3. LOADS COMBINATION USING PROBABILITY DENSITY FUNCTIONS 


Since a sine wave considered by itself is not a random signal, it has always been unclear how to 
combine harmonic analysis results with output from a random excitation analysis to obtain a design load. 
The harmonic signal can, however, be defined as a stationary random process when combined with an 
independent Gaussian random signal since its phase relationship, 0, with the random signal is random. 12 
If the harmonic signal is defined as 


/> sin =Asin(o)f + 0) , (4) 

then the argument ( (Ot + 0) is therefore a random variable with a uniform distribution over the range -rc/2 
to n/2. The probability density function of the sine function of this distribution is 


Ay) = 



(5) 


shown graphically in figure 2. 13 



Figure 2. Probability density function f(y)= sin(.v). 
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The equivalent event technique can now be used to calculate the CDF, F(z), of functions of random 
variables. 14 For the function z = x + y, where x and v are random variables, the resultant F(z) is defined in 
terms of the probability density functions (PDFs)/ v and/ v as 


oo 

fz-y \ 

F(z)= J 

j fxi x ) dx 

— oo 

V -°° / 


f v {y) d y 


(6) 


If x represents the random load, it can be characterized by a normal distribution with mean zero and 
standard deviation equal to the rms from the random analysis. Substituting the normal PDF and the PDF 
shown in equation (5) into equation (6) results in 



where the integral over y is only evaluated from -A sjn to A sin since the function is undefined outside of that 
range. A closed-form solution for this integral cannot be obtained; however, recent improvements in the 
software code Mathematica® 4.1 have enabled it to be evaluated numerically for any value of A, z, and 
(T ran . The routine crashes if the integration is carried out over the range exactly from -A to A, so the limits 
are instead set at -A+e to A-e, where e is a very small value. The value of e was obtained by performing a 
series of runs with decreasing e until the solution converged; the final value obtained was 1x10 1 ’. 

It is also of interest to obtain the solution for the inverse problem from that described above; i.e., 
what would the design load be that would yield a specific CDF value. In particular, if the 99.865-percent 
level has traditionally been chosen as an acceptable level of probability of exceedance for designing to 
random loads, then it would be consistent to use that level to design for the combination of all dynamic 
loads. The Mathematica Runge-Kutta “FindRoot” command can be implemented to calculate this load (z) 
given the probability F(z), A sin , and ( 7 ran . 
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4. LOADS COMBINATION USING MONTE CARLO 


Monte Carlo simulation can also be used to obtain the CDF value of interest. The random load is 
first simulated by generating a vector of points, { r } , that fall within a Gaussian (or normal) distribution of 
mean zero and standard deviation equal to the rms value resulting from the random analysis. This can be 
expressed as {r}~N(0.0,cr ran ). For illustrative purposes, a simulated history for such a loading is shown in 
figure 3. The vectors representing the results from the harmonic analysis are then generated. The method 
can account for analysis results from multiple independent harmonic excitations (y /=1, n. Since the 
frequencies of each sine wave are different and the relative phasing is unknown, each sine wave should be 
considered independently. 15 Examples of sine loading time histories are shown in figure 4. This method is 
only applicable for independent harmonic sources, such as excitations from different turbopumps. The 
loads coming from multiples of the primary shaft rotational speed for a single turbopump would be 
correlated and, therefore, not independent. 




Figure 4. Sample time history of sine loading. 


For the independent case, a vector for each wave of the same length as { r } is therefore created using the 
PDF defined previously, resulting in 


{y}, = A*in(2*{x},.), <*> 

where { x } is a random vector uniformly distributed between zero and 1. The resulting vector, { y } , will 
contain points that have the correct probability distribution of a sine wave, but are related randomly in 
phase (unlike the sequential time history shown in fig. 4). This procedure allows the frequencies of differ- 
ent harmonic signals to be rigorously removed from the calculation. The summation of all the harmonic 
signals and the random signal is performed for each sample to create a new vector, { z } , which is equal to a 
sample set of total responses of the structure to the multiple simultaneous excitations. 


{z} = {r} + {y},+{y} 2 + ... +{y} /f ■ 

A sample of such a loading is shown in figure 5, which is formed by adding the curves in figures 3 and 4. 
The CDF of { z } is then generated for the sample case, as shown in figure 6. A search procedure is 
performed to obtain the value from this CDF that is greater than 99.865 percent of the responses. 
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Figure 5. Sum of sine and random signals. 



Figure 6. CDF of combined loadings. 
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This procedure was implemented in a FORTRAN program running on an SGI workstation. 16 It was 
applied successfully in the MC-1 engine program in 1998 to determine design loads used to evaluate 
stresses in interface brackets, fasteners, lugs, and other structural members. Since a large number of com- 
binations were needed for each of a large number of load cycles, the speed of the program was increased by 
using efficient numerical routines found in standard FORTRAN libraries, such as IMSL, and by automat- 
ing the process for multiple combinations. Using these methods, 100,000 samples for each of 792 separate 
sine-random load combinations were calculated in a total of 1 ,573 CPU s (26 min) on a Cray T90. For more 
general use, a generic, user-friendly computer program implementing this method was also written and 
incorporated into the NASA/Techtrac national software library. 
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5. ACCURACY AND SPEED OF MONTE CARLO 


The accuracy of the Monte Carlo simulation is dependent on the number of samples selected. A 
convergence study performed here for a single harmonic wave with a random signal indicates that the 
values converge for =400,000 samples. There are several analytical methods for determining accuracy. 
Applying the method described by Ang and Tang, 17 the variation in the probability of failure estimate, p f , 
for a given number of samples, n, with a specified confidence is 


Pf 



( 10 ) 


Charts for various confidence levels are used to determine k^. For a 95-percent confidence; i.e., 95/100 
simulations would yield a probability estimate within the tolerance band, & c ^2 = l-96. The accuracy using 
400,000 samples is therefore 0.99865 ±0.0001 14 (99.854-99.876 percent). If 200,000 samples are used 
instead (doubling the speed), the error would be 0.99865+0.000161, showing only a small loss of 
accuracy. 
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6. IMPLEMENTATION IN MICROSOFT EXCEL 


To facilitate ease of use, both the PDF inverse-integration technique, discussed in Section 3, and 
the Monte Carlo technique were implemented as Visual Basic macros launched from Microsoft Excel. A 
complete listing of the macros and associated batch runstreams are presented in appendices A and B. The 
listings are also available in electronic format from the authors. In each case, any number of different 
combinations of random Id values and sine amplitudes are listed in columns, and these columns, along 
with a blank column adjacent, are chosen, the macro is selected, and then launched. The inverse-integration 
technique automatically writes the input values to a text file and launches the Mathematica kernel in the 
background, which reads the data, performs the root search for the chosen probability level (hard-coded as 
99.865 percent), and writes the result into a text file that is then read into the result column in Excel. This 
process takes =4.5 s for each combination on an Intel Pentium 4 1 .8 GHz workstation with 500 Mb RAM 
running Windows® 2000. 

The Monte Carlo macro is based upon the FORTRAN code written previously, but can only 
combine a single sine load with the random load. It is run entirely within Visual Basic for Excel to enable 
maximum portability, and takes -51s per combination. 



7. RESULTS 


The integration method was used to evaluate the actual CDF values using each of the three closed- 
form methods for a range of sine and random inputs chosen from actual MC— 1 analysis results as well as 
some hypothetical cases. A representative subset of these CDF results is shown in table 1 . The combination 
load values are also compared with loads obtained for a 99.865-percent probability level using the inverse- 
integration method and the Monte Carlo method, with the inverse-integration method being chosen as the 
baseline. The Monte Carlo method shows very close agreement resulting from the large number of samples 
arrived at from the convergence study. As expected, both the standard and the 3xssMS methods always 
exceed the baseline value by varying amounts. The peak method always yields a value less than the baseline. 


Table 1 . Comparison of all methods. 
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8. EXAMINATION OF CLOSED-FORM METHODS 


If a user chooses to apply one of the simple closed-form techniques, it would be desirable to use the 
technique that yields an answer with a CDF above 99.865 percent, but not excessively so. As seen in 
table 1, for some combinations, this rationale would point to the standard method, while for others, the 
3xssMS would apply. The exact ratio defining which is better to use can be calculated. First, equate equations ( 1 ) 

and (2): 


Asin + 3c r an “ ^(^sin)” + (°Y an ) • 

Using the relationship A sin =^2 <7 sin , equation ( 1 1 ) is satisfied when 

<T sin = 1.21218<7 ran , or A sin = (1.71 429)cr ran . (12) 

Therefore, if A sjn > 1.71429 <7 ran , the standard method will yield a probability value greater than 
99.965 percent but less than the value obtained using the 3xssMS method, and vice-versa. This result is 
demonstrated in figure 7, which is a plot of the percentage overshoot versus the ratio, -v=<r ran M sin , for each 
of the above methods. 
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The functional relationships identified in figure 7 were examined for various values of the param- 
eters and were seen to be only a function of the ratio. This independent functional relationship identified 
the possibility of creating a closed-form equation that would exactly match the values obtained using the 
inverse-integration method for this specific probability level. The Mathematica 4.1 curve fit routine was 
used to examine different categories of functions for fitting the overshoot, resulting in equation (13), with 
the comparison of the function with the data shown in figure 8: 


overshoot = 


g - 1.53546 lx 


4.60691 xl0 _11 x 7 
+1. 19197 xlO -6 * 5 
+ 1. 69965 xl0~ 3 .v 3 


v 


- 1. 1 93679 xl0 _8 x 6 
-6.041 17 xlO~ 5 x 4 
- .0271317.V 2 + .241 154.v 


\ 


! 


The desired closed-form equation is therefore 


(13) 


design load = 


3 l^sin +(J 2 

1 V2 ran 

1 + overshoot 


(14) 



Figure 8. Curve fit of overshoot of 3xssMS method over CDF of 99.865 percent. 
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A comparison of the baseline values with those obtained using this new closed-form solution is shown in 
table 2. The agreement is excellent, with only one of the points showing an overshoot of greater than 
1 percent. This equation can be directly used in a spreadsheet. If the ratio cr ran to A sjn is outside of the range 
0.025 to 10, or if either component is zero, the curve fit is not valid; therefore, an Excel user-defined 
function has been written (see app. C), which will automatically apply the curve fit within the acceptable 
range or the appropriate traditional method outside of it (where they are accurate). 


Table 2. Comparison of derived closed-form equation with 99.865-percent baseline. 


Amp. 

Sine 

la 

Random 

Ratio 

° ran/^sin 

99.865% CDF 
Using Inverse 
Integration 
(Baseline) 

Closed-form 
Curve Fit 
Based on 
3xssMs 
Method 

Overshoot 

From 

Baseline 

(%) 

1 

10 

10.000 

30.0811 

30.0749 

- 0.02 

2 

10 

5.000 

30.2986 

30.2979 

0.00 

3 

10 

3.333 

30.6500 

30.6547 

0.02 

4 

10 

2.500 

31.1169 

31.1152 

1 

CD 

O 

6 

10 

1.667 

32.3175 

32.2508 

- 0.21 

8 

10 

1.250 

33.7557 

33.5979 

- 0.47 

10 

10 

1.000 

35.3340 

35.1116 

- 0.63 

20 

10 

0.500 

44.0679 

44.0599 

- 0.02 

40 

10 

0.250 

62.8036 

63.1242 

0.51 

60 

10 

0.167 

82.0509 

82.0168 

- 0.04 

80 

10 

0.125 

101.5086 

101.2031 

- 0.30 

100 

10 

0.100 

121.0833 

120.8167 

- 0.22 

200 

10 

0.050 

219.7239 

219.6545 

- 0.03 

300 

10 

0.033 

318.9018 

319.8996 

0.31 

400 

10 

0.025 

418.3045 

417.9850 

- 0.08 

50 

1 

0.020 

51.7833 

51.7960 

0.02 
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9. CONCLUSION AND AREAS OF FUTURE WORK 


The probability values associated with several widely-used methods for combining harmonic and 
random loads have been determined. In addition, two new Microsoft Excel macros have been written that 
calculate the combined load value for a specified probability level. Because of the speed of desktop 
computers, these macros can be executed easily and quickly. A closed-form solution approximating the 
value for a level of 99.865 percent has also been determined. The comparisons presented in this TP can 
assist the analyst in choosing a method for combining these loads; obtaining an accurate value can be 
critical for the design of weight-sensitive aerospace structures. A version of the macro using Monte Carlo 
simulation was successfully applied to the determination of loads for 44 separate structural components in 
the FASTRAC/MC-1 engine program for 10 or more load cycles, resulting in large weight savings. One 
avenue of future work in this area is incorporating the effect of the correlation of the sine loads for 
harmonics of a single rotating shaft. Another will be to implement the methodology discussed here to the 
extreme-value distribution. 
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APPENDIX A— MONTE CARLO TECHNIQUE 


This appendix details the Microsoft Excel macro for performing combination of random and sine 
loads using Monte Carlo technique. Place A sin in column 1 , <x ran in column 2, leave column 3 blank, select 
all three columns, and run macro “montecarlo”. If a probability level other than 99.865 percent is chosen 
for use, this value can be inserted into line 24. 


Option Base 1 
Sub montecarlo ( ) 

maxsam = 1000000 

Dim combo (1000000) As Double 

Dim pdf (10001) As Double 

Dim div(10001) As Double 

Dim cdf (10001) As Double 

Dim cum (10001) As Double 

Dim xlo As Double 

Dim xhi As Double 

Dim delta As Double 

Dim rnorm As Double 

numsam = 400000 

Pi = 4 . 0*atn (1.0) 

cur rent_se lection = Selection .Address 

For I = 1 To Range (current_selection) . Rows . Count 

a = Range (current_selection) .Cells (I, 1 ) 
sigma = Range (current_selection) . Cells (I , 2) 

For J = 1 To numsam 

rnorm = Application . WorksheetFunction.Normlnv (Rnd, 0#, sigma) 
combo (J) = a * Sin (2 * Pi * Rnd) + rnorm 

Next J 

Call Sort (numsam, combo) 
ivalue = Int (0.99865 * numsam) 

Range (current_selection) .Cells (I, 3) = combo (ivalue) 

Next I 
End Sub 

Public Function Sort(N, RA) As Integer 18 
Dim RRA As Double 
Dim L As Long 
Dim I As Long 
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Dim IR As Long 
Dim J As Long 
Sort = 0 
L = N / 2 + 1 
IR = N 
CONTI 0 : 

If (L > 1) Then 
L = L - 1 
RRA = RA ( L) 

Else 

RRA = RA ( IR) 

RA ( IR) = RA ( 1 ) 

IR = IR - 1 
If (IR = 1) Then 
RA ( 1 ) = RRA 

Sort = 1 
Exit Function 
End If 
End If 
I = L 
J = L + L 

CONT20 : 

If (J <= IR) Then 
If (J < IR) Then 

If (RA ( J) < RA ( J + 1)) Then J = J + 1 
End If 

If (RRA < RA ( J) ) Then 
RA ( I ) = RA ( J) 

I = J 
J = J + J 
Else 

J = IR + 1 
End If 

GoTo C0NT2 0 : 

End If 
RA ( I ) = RRA 

GoTo CONTIO : 

End Function 
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APPENDIX B— INVERSE-INTEGRATION TECHNIQUE 


This appendix details the Microsoft Excel macro, DOS batch runstream, and Mathematica Kemal 
text input file for performing combination of random and sine loads using inverse-integration technique. Place 
^sin * n column I , <7 ran in column 2, leave column 3 blank, select all three columns, and run macro “eq3si^”. 

Microsoft Excel Macro “eq3sig” 

Sub eq3sig() 

current_se lection = Selection. Address 
For i = 1 To Range (current_selection) . Rows . Count 
a = Range (current_selection) .Cells (i, 1 ) 
sigma = Range (current_select ion) . Cells ( i , 2) 

Open "C:\Program Files\Wolf ram Research\Mathematica\4 . 1 
\excel_to_math.txt" For Output As #1 
Open file for output. 

Print #1, a 
Print #1, sigma 
Close #1 

shell_sync ("C:\Program Files\Wolf ram Research\Mathemat ica\4 . 1 
\mat_autorun . bat" ) 

Open "C:\Program Files\Wol f ram Research\Mathematica\4 . 1 
\math_to_excel.txt" For Input As #3 
Input #3, y 

Range (current_selection) . Cells (i , 3) = y 
Close #3 

Kill "C:\Program Files\Wolf ram Research\Mathematica\4 . 1 
\math_to_excel . txt" 

Kill "C:\Program Files\Wolf ram Research\Mathemat ica\4 . 1 
\excel_to_math . txt" 

Next i 
End Sub 

Public Function shell_sync (cmd As String) As Integer 
'shell function runs application asynchronously 

'for short duration apps where assurance of completion is desired this 
function enters a do loop activating the application until an error 
occurs signifying the app is no longer there and therefore complete 
ireturn = Shell (cmd, vbMinimizedFocus ) 

If (ireturn <> 0) Then 
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' cmd successfully initiated 
'set error trap and activate app in loop 
On Error Resume Next 
Do While Err. Number = 0 
'idle time loop 

For i = 1 To 1000 
Next i 

'check for system events that require action 
nfrm = DoEventsO 
'attempt to activate app 
AppActivate ireturn 
Loop 

'turn off error trap 
On Error GoTo 0 
'shell-sync successful 
shell_sync = 0 
Else 

'cmd initiation unsuccessful 
'shell-sync unsuccessful 
shell_sync = 1 
End If 

End Function 

DOS batch runstream “mat_autorun.bat”. Place file in directory with Mathematica Kernel executable. 
C; 

{-.hd.3.1" C:\Program Files\Wolf ram Research\Mathemat ica\4 . 1 
Math crandsin calc.txt >math_output . txt 


Mathematica Kernel text input file “randsin_calc.txt”. Place file in directory with Mathematica Kernel 
executable”. If a probability level other than 99.865 percent is chosen for use, this value can be inserted 

into line 6. 

AppendTo [$Echo, "stdout"] 

Off [General : : spelll] 

v= Import ["excel_to_math.txt" , "List"] ; 
a=v [ [ 1 ] ] ; s igma=v [ [ 2 ] ] ; 

\ ! \ (max = 3*sigma\ + a; eps=l . 0* 10 A - 11 ; min = ,85*max;\n 

FindRoot [ .99865 == \(l\/\(a\ * 

sigma \ \@\(2 \ [Pi] \ A 3\) \) \) \ (\ [Integral] \_\ ( -\ ( (a - 

\eps) \) \) \%\ ( (a - eps) \) \ { (\ [Integral] \_\ (\ ( -6\) *sigma\) \%\ (z 

\ y \ ) \ ( \ [Exponent ialE] \ A \(\(-\(( X \/ sigma) \)\ 2\)\/2\)\) 

\ [Dif f erentialD] x) \ ) \ 
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\(1\/\@\(1 - \ ( (y\/a) \) \ A 2\) \) \ [DifferentialD]y\) , {z, {min, max} } ] \ ) 

equiv3sigma=z/ . % 

Export [ "math_to_excel . txt" , equiv3sigma, "List "] ; 

Quit [ ] ; 
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APPENDIX C— APPLYING CURVE FIT OF THE PROBABILITY DENSITY FUNCTION 


This appendix details the Microsoft Excel user-defined function applying curve fit of PDF Integra- 
tion. Create macro as shown below, then in spreadsheet, type function “curvefit(a,s) with A sin 
as the first argument and <T ran as the second argument. 

Function curvefit(a, s) 

If a = 0 Then 

combo = 3 * s 
Elself s = 0 Then 
combo = a 

Elself (s >= 0.02 * a) Then 
x = s / a 

combo = ( 3 * Sqr (a2/2+s2))/_ ^ 

(1 + Exp ( -1 . 53546*x) * ( (4 . 60691E-ll/x A 7) - (1 . 19367935364834E-08/x 6)+ _ 

(0.00000119197 /x*5) - (0.0000604117 /x 4) + (0.00169965 /x 3) - _ 

(0.0271317 /x a 2) + ( 0 . 24 1154/x) ) ) 

Else 

combo = a + ( 3 * s ) 

End If 

If (s > (10 * a) ) Then 

combo = 3 * Sqr (a 2 / 2 + s 2) 

End If 

curvefit = combo 
End Function 
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